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Abstract 

We present the first search for the decay — t- we"*"!/ to test the four-quark content of tlic 
and the u-cj) mixing model for this decay. We use 586 pb~^ of e+e~ colhsion data collected at a 
center-of-mass energy of 4170 MeV. We find no evidence of a signal, and set an upper limit on the 
branching fraction of B{Df ue^u) <0.20% at the 90% confidence level. 



* Now at: Pacific Northwest National Laboratory, Richland, WA 99352 
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I. INTRODUCTION 



Multiple observations of exotic cliarmonium states [1 — 4] have been widely interpreted 
as four-quark states [5 — 13]. In this analysis we probe the four-quark content of the 
by searching for the decay — )■ ue'^u (charge conjugate states are implied throughout the 
article). Assuming that the u is a pure two-quark state, its valence quarks are distinct from 
those of the D^, and the decay can proceed through the diagram of Fig. 1, where the (uu) or 
(dd) come from within the D^. The initial valence quarks annihilate while a lepton pair is 
produced. Neither Cabibbo-favored, nor Cabibbo-suppressed decays can contribute to this 
final state. The study of this specific process was first suggested in Ref. [13], and Ref. 15 
estimates the theoretical branching fraction for the analogous decay J/ip i~^i>. 




UO 



FIG. 1. Feynman diagram representing the four-quark semileptonic decay Df — )• uje'^v. 



Recent work by Gronau and Rosner pi)| concludes that any value of the branching fraction 
for Z)+ — 7- ue^v exceeding 2 x 10~^ is unlikely to be explainable by u-cj) mixing and would 
provide evidence for non-perturbative effects known as "weak annihilation" (see Ref. [ii 



for references). An estimate based on comparing hadronic and semileptonic processes gives 
a branching fraction of (0.13 ± 0.05)%. This is the first search for this decay mode. 

We search for a positron candidate and an w — 7r+'7r~7r° candidate, which is the dominant 
decay mode with a branching fraction of 89.2% [17|. Cabibbo-favored decays exist in the 
same final state, — )■ r^e+z/ and — )■ (j)e'^i', with B{r] — )■ 7r"'"7r~7r°) =22.73% and and 



S(0 — > 7r+7r^7r°) =15.32% 17|. They can play the role of control samples, which are used 
directly in the analysis in a variety of ways. For example, the effect of certain selection 
requirements can be readily estimated from any change in the rj and populations. The 
two control samples are also well measured using the independent final states 77 — )■ 77 and 
— )• K^K~ . Therefore, this analysis has good statistical sensitivity, and unusually strong 
control samples using CLEO-c data directly. 

The remainder of this article is organized as follows. In Sec. II, the detector, data, and 
Monte Carlo (MC) samples are described. The data analysis method is described in Sec. III. 
The fitting procedure is described in Sec. IV. The determination of the branching fraction 
is discussed in Sec. V. 
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II. DETECTOR, DATA AND MC SAMPLES 



The data used in the analysis are e^e^ colhsions at a center-of-mass energy ^/s =4170 
MeV. The data are collected by the CLEO-c detector and correspond to an integrated 
luminosity of 586 pb~^, or 0.6 x 10® D^D~ inclusive pairs. 

The CLEO-c detector is optimized for physics in the charmonium region and is described 
in detail in Ref. [lH. The tracking system consists of a central, low mass drift chamber, 
which wraps directly around the beam pipe, and a main drift chamber inside a solenoidal 
magnetic field. The particle identification system combines the track information from the 
gas chambers (dE/dx) and the associated ring- imaging Cherenkov detector (RICH) data. 
The electromagnetic calorimeter consists of Csl crystals, arranged in cylindrical fashion 
around the drift chamber, to make a barrel at angles given by | cos^| < 0.7, with 6 being 
the angle measured from the interaction point (IP) with respect to the beam axis. The ends 
of the cylinder are also instrumented with Csl crystals and are referred to as the end-cap 
regions. 

MC simulations of the known physics processes j20l] and of the CLEO-c detector are 
used to estimate backgrounds and calculate signal resolution and efficiencies. The known 
charm physics processes are included in the cc MC simulation. All types of charm back- 
grounds, dominant in this analysis, are simulated to 20 times the statistics in the data, while 
the continuum {u, d, s) backgrounds are simulated to 6.6 times the statistics in the data. In 
the following, where MC results are presented, we multiply the continuum MC sample by 3 
to obtain a consistent x20 normalization. By convention, cc MC refers to the charm part 
of the MC, continuum MC is the non-charm part, and MC is the weighted sum of the two. 

The signal MC sample consists of 8x10^ — )■ ue^u events generated according to 
phase space distribution. The same number of events was generated for the — rje^u 
and — )■ (pe'^u control samples. 

Further samples for dominant sources of backgrounds were generated. A sample of 2x 10^ 
— !■ (jj7r+7r° events was separately generated, corresponding to about 23.5 times the num- 
ber of such expected decays in the data. In addition, 5000 events (96 times the data) were 
generated for the decay chain — )■ rj'e^v followed by rj' — )■ 077. 



III. DATA SELECTION 

About 95% of the e"'"e~ — )■ DfX events are formed through the following exclusive 
reactions 18 

e+e" ^ D^D^- and e+e" ^ A"^s^- (1) 



Equal amounts of positive and negative D* states are produced, and about 95% of them 
decay to Ds7- The analysis described here selects exclusively both final states in Eq. ([1]). 
In the following, we retain both events where the 7 is associated with the positive side, the 
signal side due to the convention established in Sec. I, and where the 7 is associated with 
the negative side, or tag side. Several kinematic constraints are available, but only those 
that select both reactions in Eq. ([T]) are used. 

First we search for an exclusively reconstructed hadronic D~ candidate, the tag, and a 
photon candidate. Requiring both of these objects, with three kinematic selections applied, 
strongly suppresses the backgrounds. The photon candidate, and all tracks and showers 
which are daughters of the tag candidate, are not used when searching the rest of the 
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event, which, due to the exclusive nature of the analysis, must be the signal candidate. 
The following are required to form a signal candidate: a positron of opposite charge to 
the tag, precisely three charged tracks, a 7r° candidate, and net event charge equal to zero. 
Furthermore, the missing energy and momentum are required to be in a relation consistent 
with the presence of a nearly massless neutrino. Extra showers in the event are ignored. 
The selection is described in detail in the remainder of this section. 



A. Charged and neutral particles selection 

The criteria for a good track include the requirements that its minimum distance to the 
IP does not exceed 5 mm in the plane perpendicular to the beam axis and 5 cm along the 
beam axis. Phase space is limited to | cos^| < 0.93, with 9 being the angle with the beam 
axis, and momentum 0.05 GeV < p < 2 GeV. Good tracks are then selected as pion or kaon 
using dE/dx and RICH data according to the algorithms described in detail in Ref. (22| . 

Photon candidates are contiguous groups of crystals recording significant energy deposi- 
tion. They are required to be unmatched to tracks or noisy crystals, and to have a transverse 
profile consistent with expectations from an electromagnetic shower. The minimum cluster 
energy is 30 MeV. 

Positrons are selected by requiring a good track, and the combined positron probability 
of the particle ID system and E/p (the ratio between the shower energy associated with the 
candidate electron and its track momentum) to be greater than 0.8. Positron phase space 
requirements are | cos^| < 0.9 and p > 0.2 GeV. 

7r° and 1] candidates are selected by requiring a photon candidate in the barrel or end- 
cap, with I cos^l < 0.93. This photon is combined with a second photon, which must not be 
associated with a noisy crystal, and the invariant mass of the two photons must be within 
3 standard deviations of the nominal 7r° and t] masses [l7| . 

candidates are selected by requiring two oppositely charged tracks. If they are assigned 



the nominal vr"*" mass 17|], their invariant mass must be within 12 MeV of the nominal 



mass [17[. A common vertex is calculated, and it is required to be radially displaced from 
the IP by at least 3 standard deviations. 



TABLE I. Definitions of Mtag signal and sideband definitions for the tag modes. 



Mode Signal region (GeV) Low sideband (GeV) High sideband (GeV) 



KlK- 


[L954, 


1.983] 


[1.910, 


1.939] 


[1.998, 2.026] 




[L954, 


1.982] 


[1.911, 


1.940] 


[1.996, 2.025] 


K*-K*° 


[L953, 


1.983] 


[1.909, 


1.938] 


[1.997, 2.027] 




[L955, 


1.982] 


[1.913, 


1.941] 


[1.996, 2.024] 




[L940, 


2.001] 


[1.892, 


1.922] 


[2.019, 2.050] 


r]p- 


[L940, 


1.998] 


[1.885, 


1.914] 


[2.021, 2.050] 




[L944, 


1.992] 


[1.885, 


1.933] 


[2.004, 2.052] 




[L944, 


1.992] 


[1.886, 


1.930] 


[2.002, 2.047] 
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B. Tag-candidate selection 



Eight tag decay modes are used and listed in Table [U using the particle candidates 
selected according to Sec. III. A. In addition, there are several mode-specific criteria. For 
— i- K^K^TT^ and tt+tt^tt^, the pion momenta are required to be greater than 



0.1 GeV. For — )■ K* K , only the {K^tx ){K^tx ) channel is considered. The K* and 



K candidate masses are required to be within 100 MeV of the nominal value [17|. For the 
D~ — i- rjp^i^p" — )■ 7r~'7r°), the p~ mass must be within 150 MeV of the nominal value 17 . 
For -Do—?- Ti^vlirl — > p^\p^ — )■ 7r"'"7r~, the Vj mass must be within 20 MeV of the nominal 
value 13]. Furthermore, the [71-1]') helicity angle (defined as the angle 6h, in the rest frame 
of the p, between the momentum of the vr" and the momentum of the D~) is required to 
satisfy | cos 9h\ < 0.8. 

The four-momentum of a tag candidate is defined by (-Stag, Ptag), with the tag mass 
defined by Mtag. The selection further makes use of the recoil mass Mrec, defined as 



M,ec = V(^b - ^tag)2 - (Pb - Ptag)2. (2) 

Here, (i?b, Pb) is the four-momentum of the colliding beams. The Mrec distribution will 
peak only for those events where the photon is associated with the signal side, but even 
when the photon is associated with the tag side, Mrec is kinematically constrained so that 
I Mrec ~ M*| < 55 MeV, where M* is the nominal D* mass fv?\. Only candidates passing 
this selection are retained. 

The main tag selection is obtained from a 2-D fit described below. Because of the com- 
plexity of the fit, the Mtag projection is fitted first, and the fit results used to constrain 
some of the final 2-D fit nuisance parameters. The Mtag projection is also best suited for 
side-band background subtraction. 

The Mtag distribution is fitted with a double Gaussian function, G2, multiplied by the 
fitted number of events, A^, and a first degree polynomial, Ai, to describe signal and back- 
ground 

/(Mtag) = ArG'2(Mtag) + Ai(Mtag). (3) 

G2 is a probability distribution composed of two Gaussian functions G{x; a, fi), of unit area, 
peaking at fi and with width equal to a. Here, the peak is set at Mds, which is the nominal 
nominal mass [13], and the two Gaussians have fractional probabilities f\ and (1 — /i) 

G'2(Mtag) = /iG(Mtag; 0-1, Mds) + (1 - /i)G(Mtag; cxa, Mds). (4) 

The quantities cxi and a2 are fixed to the value obtained from the fit to the signal MC 
data. Having obtained /i from the fit, we construct the variable (7^2 = ficf + (1 — /i)cr|. 
The signal regions are required to be within 2.5cti2 from the peak position for each mode 
except for the [rip) mode where it is selected within 2cri2. 

The Mtag data fit results are listed in Table [Tll The rest of the peak fit parameters 
are listed in Table IIIII The Mtag distributions for the 8 modes are shown in Fig. [2l The 
sidebands in the Mtag distribution are listed, for each mode, in Table [H 

Having determined the fit parameters for the Mtag distributions, a second kinematic 
constraint can be imposed using the MM*^ variable defined as 

MM*2 = (Eb - Etag - E^f - (pb - Ptag - P7)', (5) 
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where (i?^,p^) is the photon four- momentum. If the final state is given by Eq. ([T]), MM*^ 
should peak at M^^. The MM*^ mass selection criteria are found by a two-dimensional 
(2-D) binned likelihood fit in the (MM*2,Mtag) space. Each variable is also kinematically 
fitted, so that Mtag is the value obtained by constraining MM*^ to its nominal value, and 
vice versa. This procedure improves the signal and also minimizes any correlation between 
the two variables. 

TABLE II. Number of Dj tag candidates A'^data for each mode in the signal region and sidebands 
for the one-dimensional fit to Mtag. 



Modes N^ata Low Sideband High Sideband 





5828 ± 92 


1231 


958 


K+K-JT- 


25990 ± 285 


22385 


19452 


K*-K*° 


2891 ± 100 


2783 


2647 




8152 ± 369 


56530 


43475 




3635 ± 160 


5727 


3379 


rip- 


6877 ± 330 


26879 


14658 


7r~r/'(r/7r^7r~ 


) 2344 ± 70 


1040 


572 


7r-r/'(p7) 


4451 ± 337 


42412 


25476 



TABLE III. Signal peak parameters for each tag mode derived from cc MC simulation. 



Mode /i fJi (MeV) tJa (MeV) 

K^K- 0.471 405 7!00 

K+K-TT- 0.725 3.74 8.92 

K^-K*^ 0.771 3.43 10.65 

vr+TT-vr- 0.899 4.84 9.88 

riTT- 0.650 9.85 15.56 

rip- 0.574 10.8 18.3 

TT-r]'{r]7r+TT-) 0.590 5.71 13.34 

TT-'n'{p-f) - 9.60 



The 2-D fit is done for each mode separately, and its purpose is to extract the final number 
of tags for each mode, Ni, while building on the information obtained in the one-dimensional 
Mtag fit. The fitting function is 

f{M,,^,MM*^) = Ar,G2(Mtag)C(MM*2) + G'2(Mtag)A5(MM*2)+Ai(Mtag)A5(MM*2). (6) 

For each mode, the signal is described by the product of a double Gaussian in the Mtag 
projection, defined in Eq. (jlj) and Crystal Ball function jlsj, defined as C in the equations, in 
the MM*"^ projection, respectively. One of the background components is the combination 
of a real tag with a random 7. This type of background {BGi in Fig. 3 below, and the 
second term in Eq. ([6])) is described by the same double Gaussian G2(M;ag) and a 5th degree 
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FIG. 2. Distribution of Mtag — Mds of candidates, for the different tags: (a) K^K ir ; (b) 



K^K~; (c) r?vr~; (d) tt^t]' {rj-K^-K^); (e) vr+vr^vr"; (f) K*~K* ; (g) r}p~ , and (h) it~ rj' (p'j) . The 
fitted background Ai, described in the text, is indicated by the dashed-dotted slope. The signal 
mass region is indicated by the vertical dotted lines. 



polynomial A^lMM*"^). The other background {BG2 in Fig. 3 below, and the third term in 
Eq. ([6])) is due to fake tags. The PDF here is the product of a 1st order polynomial (Mtag) 
and a 5th order polynomial (MM*^). To simplify the fit, the Mtag projections are fitted 
using the signal function obtained in the 1-D Mtag fit, but the background parameters are 
varied. 

TABLE IV. MM*^ selection range and the number of tags in each mode obtained from the two- 
dimensional fit. Total number of tags, -/Vtag is also given. The quoted error is statistical only. 



Modes 


Lower hmit ( GeV^) Upper limit ( GeV^) A''i(data) 


K^K- 


3.7876 


3.9539 


3442 ± 138 


K+K-TT- 


3.7939 


3.9510 


15647 ± 271 


K*-K*^ 


3.7505 


3.9847 


1707 ± 94 




3.7701 


3.9633 


4595 ± 298 




3.7662 


3.9798 


2355 ± 187 


VP" 


3.7698 


3.9632 


3606 ± 640 




) 3.7409 


3.9888 


1716 ± 142 




3.7875 


3.9601 


3373 ± 240 


iVtag 






36441 ± 852 



The MM*^ distributions are shown in Fig. [31 The MM*^ signal regions for each mode are 
chosen so as to have 95% signal efficiency. The MM*^ selection is summarised in Table IIVI 
Table IIVI also lists the final number of tags obtained in each tag mode, A^j, as well as the 
total number of tags, A^tag, used to extract the final result. 
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(GeV)2 

FIG. 3. MM*^ distributions for the 8 tag modes. Dash-dotted Hnes are the BG2 background 
described in the text. Dashed hnes are total background BGi + BG2- (a) ; (b) Kq^K~; 

(c) rjT[~; (d) TT~r]' {r]7r~^TT~); (e) vr+Tr^vr^; (f) K*~TC*^; (g) rjp^ , and (h) 7r^r/'(p7). 



C. Signal selection 

The signal is selected by requiring one positron candidate, of charge opposite to the tag 
charge, two charged pion candidates, of opposite charge, no extra good tracks, and a good 
7r°, all selected exclusively of the objects used in the tag. The selection requires a specific 
number of tracks, and multiple candidates can arise only due to multiple 7r° candidates. In 
case of multiple candidates, the 7r° is selected as follows. Given the photon-photon mass 
My^, and the calculated mass error cr^^, the one with the lowest = [(^77 ~ M^^o) / a-y-y]"^ is 
chosen. Additional candidate photons are ignored. 

The positron, charged pions, and 7r° are added together to form the four- vector (£'s,Ps). 
The measured neutrino candidate mass squared, MM^, is defined as 

= (Eb - Etag - ^7 - ^s)' - (Pb - Ptag - P7 - Ps)'- 

The MM^ distributions, with Mtag sideband subtraction, of the two control samples rje'^u 
and (pe'^u are shown in Fig. |H Based on the shape of MM"^, events with —0.05 GeV^ < 
MM^ < 0.05 GeV^ are selected for the final analysis. 

The mass of the tt'^h'tt^ combination M3 was not used in the candidate selection, and 
provides the spectrum that is fitted to extract the final result. In Fig. the M3 spectrum 
is presented, including Mtag sideband contributions. Two peaks are clearly present, at the 
rj and (p masses, with no sign of a signal in the u mass region. 

The M3 peaks in the signal MC samples are fitted to a Breit-Wigner shape (indicated as 
BW in the equations), convoluted with a double Gaussian, 

s{x) = K J BW{xi)G2{x - xi)dxi, (7) 

where K is a. normalization constant to give s{x) a unit area. Table |V] lists the fit results 
for each of the signal MC samples generated for this analysis. 
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FIG. 4. Control sample MM"^ distributions, (a) solid, D'^ — )• rje^v distribution after Mtag sideband 
subtraction, and Mtag sideband distribution (dotted), (b) solid, Z)+ — )• -qe^v distribution after Mtag 
sideband subtraction, and Mtag sideband distribution (dotted). 



TABLE V. M3 signal peak parameters evaluated from signal MC sample. All quantities are defined 
in the text. 



Decay 


fi 


C7i (MeV) 


a2 (MeV) R.M.S. (MeV) 




0.8844 


3.165 


19.85 7.37 




0.8783 


5.500 


22.52 9.40 




0.8361 


5.940 


19.83 9.73 



The reconstruction efficiency e for the ue^u final state is computed by applying the same 
requirements to the signal MC events, but correcting for the number of tags found in the 
data, 

e = ^SiV,e„ (8) 

-'^tag 

ej being the signal MC efficiency for mode i. The result is e = (5.11 ±0.15)%, with the error 
due to MC statistics. 



IV. FINAL FIT 



Figure [6] shows only the M3 region used in the fit, which contains A^^obs = 18 events. The 
Ajv/3 = 250 MeV mass window is centered at the nominal u mass 17|]. In the Fig. 6(a), the 
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FIG. 5. M3 distribution. Solid: signal selection, after Mtag sideband subtraction. Dotted: Mtag 
sideband contribution. The arrow shows the location of the oj nominal mass, Ref. [17]. 



data distribution and the fit to the data (described below) are shown. In Fig. 6(b), the MC 
distribution is shown. 

Three potential sources of background are considered: non-Dg backgrounds, back- 
grounds where there are non-resonant final states (which have not yet been observed, and 
are not present in the MC simulation), and backgrounds where there is a true u. Mtag 
sideband subtraction only subtracts the first source. A direct fit of a signal and a back- 
ground component subtracts the first two. The third source of background is subtracted via 
MC simulation, and is discussed below. The signal yield is determined by a one parameter 



unbinned likelihood fit [17[. The free parameter is the total number of signal events S. 
The background level is constrained by the normalization of the probability. S is multiplied 
by a function of unit area s{x), Eq. (|7]). The final expression of the unbinned likelihood, 
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FIG. 6. AI3 distribution in the 250 MeV wide region centered at the nominal uj mass. No Mtag 
sideband subtraction was used, (a) comparison of data (points) and best fit according to Eq. (9) 
(hue), (b) The MC sample, normalized to 20 times the data statistics. Solid: non-resonant 
backgrounds. Empty: resonant backgrounds. The arrow shows the location of the cj nominal 
mass. 



Lu = HjPj, is obtained from the probabilities 

P{M3i\S) = Pi = {S/N,^,)s{M,i) + (1 - S/N,^,)/Am3, (9) 

which correspond to a signal S, distributed according to s^M^), plus a fiat background. 

Figure shows the likelihoods obtained for data, without any peaking background sub- 
traction, in the S > region. The 90% confidence level (CL.) is calculated using only the 
S > portion of the likelihood. The statistical only upper limit on S at the 90% CL. is 
5*90 = 3.78 events. 

Unbinned likelihood fits, in one dimension, can be tested for goodness of fit using the 
Cramer- Von Mises test [24!], where the goodness of fit parameter is 

G= [F{Ms)-FNiMs)]^dF{Ms). (10) 

J A/™'" 

The integral limits are the limits of the fit interval. FlM^) is the integrated probability 
function for best-fit parameters, 

FiMs)= , P{M!„S^,,)dM;,. (11) 

Here Smax = 0, so that F is in fact a straight line. One has F = (M3 — M™™)/Am3 and 
dF{M^) = dMs/Ajvfa- F^- is a step function such that Fj^{M^) = N/Nohs, where is the 
rank of the largest event mass which is less than M3. The two functions are shown in Fig. [HI 
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FIG. 7. Statistical only likelihood, normalized to unit area over the positive signal region, for the 
observed number of signal events in data. The 90% C.L. was computed using only the positive 
signal region. 



A toy MC program was run to generate an ensemble of 10^ unbiased experiments. Figure [8] 
shows the distribution of G for the ensemble, also shown is the value of G obtained in the 
fit to data. Only 13.1% of the fits to the generated experiments are better than that to 
the data. The toy MC program also made it easy to apply the Kolmogorov-Smirnov (KS) 
test, which simply computes the maximal difference between F and Fj^. Only 16.0% of 
the unbiased experiments produced a better KS test than the data. The fit to the data is 
excellent. 



V. DETERMINATION OF BRANCHING FRACTION AND SYSTEMATIC ER- 
RORS 



The statistical upper limit on the number of events is translated into a statistical only 
limit on the branching fraction Bqq according to the following equation 



90 



'90 



(12) 



tag 



There are three quantities on the right hand side of Eq. flT2|) . with central values Si 



90 



3.78, e = 0.0511 and iVtag = 36441, yielding B^q = 0.203%, which is a purely statistical 
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3540311-008 




G(10-2) 



FIG. 8. Cramer- Von Mises test of goodness of fit, for the final fit of this analysis, (a) Comparison 
of the integrated probability distribution F (dashed) and the step function described in the 
text (solid), (b) Comparison of the G obtained in this analysis, with a distribution obtained from 
10^ unbiased toy MC fits. 



limit. A^tag has a statistical error of its own, and each of the three quantities in Eq. (fT2l) has 
systematic errors which are discussed below. 



A. Systematic errors 

Table IVll first row, contains the relevant parameters of the unbinned likelihood (Fig. [7j), 
in the form /i ± cr. /i is the S value for which is maximal, if one allows also S" < values. 
It describes the form of the likelihood, but is not used in the determination of the final 
result. 

Systematic errors to Sgo are also listed in Table IVTl The error associated with the assumed 
mass and width of uj is estimated by varying the central values by the uncertainties given 
in Ref. [T^. 

The greatest source of S'go systematic errors is related to irreducible backgrounds. These 
also shift the location of the likelihood peak to lower values. Fig. |6] shows the background 
distribution by physical source. All but one of the true u are due to the decay chain 
D+ ^ ie+v with [B = (1.12 ±0.35)%], followed by rj' u'y, with [B = (3.02 ±0.33)%] [TJ. 
A dedicated MC simulation for this channel generated 5000 events of which 51 passed all 
selections. This corresponds to an irreducible background of (0.53 ±0.19) events. Note that 
the largest source of error is the semileptonic B error from Ref. [l^. Therefore, this source 
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of systematics can not be significantly improved with increased simulation statistics. 

A second irreducible background comes from Dg — >■ uX events. Zero events are found in 
the MC events from the direct decay — )■ wtt^. The cc MC significantly underestimates 
the {Ds — )■ uX) yield, which is 0.6% in the cc MC but 6.1% in data [25|. The decay 
—7- UTi^ is in the cc MC, but no other u)n{7c) decays. A dedicated MC for — )■ utt^tt^, 
which was assumed to saturate the 5.5% difference, was run, and zero events were found. 
The probability for n background events, given zero MC candidates, is exponential in shape. 
The systematic error from this source can be represented as (0.02 it 0.02) in Table IVll 

There was one more true u event which is in the continuum MC sample, corresponding 
to one more irreducible background of 0.15 ± 0.15 events. 

TABLE VI. Summary of statistical and systematic errors of SgQ. The first block is the statistical 
error from the experiment. The second block are errors associated with the quantity Sgo of Eq. 
(12) and consists of uncertainties due to the mass and width of the omega, and errors in estimating 
three irreducible backgrounds as described in the text. The last block are percentage systematic 
errors associated with A'tag or e. 



Type 


Cent. val. (evts.) 


a (evts.) 


Data fit 


-0.25 


2.21 


CO mass 




0.04 


CO width 




0.006 


— ^ r/'e+i^ 


-0.53 


0.19 


Dt ^ coX 


-0.02 


0.02 


Continuum 


-0.15 


0.15 


Type 


Cent. val.(%) 


a{%) 


A'tag Stat. 




2.3 


Nti,g syst. 




2.0 


MC statistics 




2.7 


B{co Svr) 




0.8 


Tracking 




0.9 


vfO eff. 




1.0 


vr*^ selection variation 




0.5 


Positron eff. 




0.6 


MC form factor 




0.5 


Extra track selection 




0.04 


Particle ID 




0.1 



A^tag was obtained through a fit, with a statistical error of 2.3%. Systematic errors can 
enter the analysis only through the bias in the choice of fitting function. This can be 
quantified by varying the fitting function. For each tag mode, the fitting function for the 
signal was changed, first term of Eq. ([6]). The Crystal Ball function was varied in two 
ways, by keeping the n parameter fixed to its MC fitted values and by changing the (n, a) 
parameters by one a in a mode specific way. The background was also varied. Instead of a 
fifth degree polynomial, the data were fitted with a fourth and a sixth degree polynomial. 
The background was also changed by fixing the amount of BGi background (described in 
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Sec. III.B) to one standard deviation above or below its central value. Variations of A^tag 
due to changes in the fitting function were as low as —1.8% and as high as +1.5%. The 
assigned systematic A^tag error is 2.0%. 

Correlations may affect the fit of the (Mtag, MM*^) peak, because the fit assumes the two 
variables are not correlated. To study this, we have computed the (Mtag, MM*^) correlation 
coefficient in the signal MC sample, by calculating the correlation coefficient in each tag 
mode, and then reweighting for the observed number of events in each mode. The result is 
p = (—1.6 ±0.9)%. The fit error due to remnant correlations is of order and is neglected. 

Finally, there are the systematic uncertainties on the efficiency to be considered. The cu 



branching fraction uncertainty is 0.8% [17[. The tracking efficiency error is a 0.3% Gaussian 



systematic error per signal track, to be added linearly, totaling 0.9% per event |26 . 

The 7r° reconstruction efficiency error is 1% [23], but depends on the exact selection 
criteria. To estimate the size of the systematics induced by changing selection criteria, the 
signal MC sample with and without the energy and angular criteria which were used to 
select photon candidates. There were 41269 reconstructed events with the criteria, and 
41868 without the criteria, a difference of 1.5%. There were 101 events instead of 99 in the 
combined rj and peaks, a difference of 2%. We assumed a further 0.5% systematic error, 
listed in Table |VT1 



The positron reconstruction efficiency is evaluated in a manner similar to Ref. [27 
Positron efficiencies have been investigated by the Collaboration using a variety of well- 
known kinematically constrained QED processes. The experimentally measured corrections 
are convoluted with the positron momentum distribution to obtain the efficiency uncertainty, 
which is 0.6%. The effects of the extra track cuts and of particle ID cuts can be estimated 
from the MC sample, by varying or eliminating the cuts. We find errors of 0.04% and 0.1% 
for the extra track cut and particle ID cuts respectively. 

Finally, a different form factor will change the efficiency, mostly because events with low 
(the positron-neutrino mass) produce lower energy positrons. The signal MC produce 
phase-space distributed events, and therefore a constant form factor. To evaluate this source 
of systematics, the form factor was varied by ±20%, by reweighting the signal MC events 
according to the weights 

2 0.2{qf- < >) 
w±i{qi ) = 1 ± -^-^ , 

with < q^ > being the mean in the signal MC sample. The new efficiencies are 5.08% 
and 5.14% respectively, to be compared to the given value of 5.11%. A systematic error of 
0.5% is assigned to this systematics. 



B. Determination of the branching fraction 

To obtain the final result, Gaussian and non-Gaussian errors are convoluted with the 
non-Gaussian signal S distribution given by the likelihood (Fig. [7]), by means of a toy 
MC program. The procedure is used for several reasons. First, the main source of error, 
the unbinned likelihood, is non-Gaussian, whereas the smaller sources of error are mostly 
Gaussian. Second, some sources of error shift the central value of the likelihood, an effect 
which can be treated exactly by shifting the likelihood on an event-by-event basis. Finally, 
the exponential and correlated nature of some of the error sources can be reproduced exactly 
by MC simulation. 
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FIG. 9. Final probability distribution for the branching fraction, after convolution of all statistical 
and systematic errors, as described in the text. 



A total of 25 X 10® toy experiments are generated, to obtain the probability distribution 
in Fig. ini The cumulative effect of the systematic errors is to increase the limit, and the 
cumulative effect of the irreducible backgrounds is to decrease the limit. Taking into ac- 
count the systematic uncertainties and the irreducible backgrounds, the upper limit on the 
branching fraction changes from 0.203% (statistical only) to 0.201%. 



VI. CONCLUSION 

We report the first measurement of an upper limit for the branching fraction B{D'^ — )■ 
ue+u). We find B{D+ ue+u) <0.20% at the 90% C.L., which does not exclude that 
expected from the model of Ref. [l6[. 
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